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Abstract — An analysis of the characteristic function of 
Gaussian quadratic forms is presented in 1 to study the 
performance of multichannel communication systems. This 
technical report reviews this analysis, obtaining alternative 
expressions to original ones in compact matrix format. 

Keywords — Quadratic forms, characteristic function, mul- 
tichannel communication systems. 



I. Introduction 

STATISTICS of complex Gaussian quadratic forms are 
a powerful mathematical tool when dealing with digi- 
tal communications problems. Currently, they are widely 
used to analyze the performance of digital modulations over 
fading channels [2]- |S]. In [1] is presented an analysis of 
one statistic of complex Gaussian quadratic forms which is 
a very interesting and useful tool to calculate error rates. 
This analysis, based on the characteristic function of Gaus- 
sian quadratic forms, provides a set of expressions which is 
applicable to a good deal of digital communications prob- 
lems. To provide a few examples, these expressions have 
been used to obtain closed-form bit error rate (BER) results 
in fading channels for multichannel binary signals [2 -[5 ], 
quadrature-amplitude modulation (QAM) [BJ, orthogonal- 
frequency-division- multiplexing (OFDM) [7], and non- 
orthogonal multi-pulse modulation (NMM) 8 . In addi- 
tion, the results presented in [I] are also used to analyze 
other statistics of complex Gaussian quadratic forms in [5] . 

In this document is performed a revision to the analysis 
of the characteristic function of complex gaussian quadratic 
forms presented in [T|. Alternative expressions in compact 
matrix format are derived and little mistakes in the original 
ones has been detected. Up until now, research based on 
the original expressions has yielded valid results because 
it dealt with cases where the original and corrected ex- 
pressions give the same results. However, future uses of 
the original expressions could lead to wrong results. After 
deriving alternative expressions, a simple example is pre- 
sented to confirm the validity of the corrected expressions. 

The remainder of this document is organized as follows: 
Section [TT] states the general problem dealt with herein; 
Section IIIII describes the alternative analysis which gives 
the new expressions and detects the mistakes in the orig- 
inal expressions; Section IIVI presents the aforementioned 
example which confirms the validity of the new ones; and 
conclusions are provided in section fVl 



II. General problem statement 

The decision variable at the detector of several communi- 
cation systems, e.g., those employing multichannel binary 
signals |TJ , can be expressed as a sum of quadratic forms 
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where Q is a Hermitian matrix and X k , Y k complex 
circularly-symmetric Gaussian random variables. The set 
of vectors {z k } k=i,...,L are mutually statistically indepen- 
dent with a different mean but identical non-singular 
covariance matrix R, respectively defined as 
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R^E{(z fc -m fc )(z fc -m fc )«}. (3) 



The constants A, B and C must be appropriately iden- 
tified with the specific parameters of the problem, so that 
D characterizes the decision variable at the output of the 
detector to further calculate the probability of error as 
Pr{Z? < 0}. This calculation leads to non-trivial results 
when Q is indefinite, i.e., when |C| — AB > 0, because 
otherwise the probability Pr{D < 0} will be either or 
1. Note that indefiniteness for 2 x 2 matrices implies non- 
singularity. 

An expression for the probability Pr{D < 0} is derived in 
PQ using the characteristic function of D, denoted as 4>d{v). 
In this derivation, mutual independence of the terms d k 
in ([I]) is used, thus, the function 4>d(v) can be expressed 
as the product of the characteristic functions <f> k of each 
summand, i.e., 
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The first two columns of tabic U lists the parameters v\ , 
V2, OL\ki 012k-, as well as other intermediate definitions and 
the expression of the probability Pr{D < 0} derived in [TJ. 
In the table, Qi and I n are the first-order Marcum Q- 
function and the modified Bessel function of the first kind, 
respectively. 

III. Alternative analysis 

An alternative formulation of the characteristic function 
4>k is used in this section to obtain new expressions for the 
probability Pr{D < 0}. 

In |10) . Turin deducts the following expression for the 
characteristic function of a quadratic form of n variables 
z w Qz 

exp (-m^R- 1 [I - (I - jvRQr 1 } m) 

^ v) = \^]vm\ ' (6) 

where Q is an n x n Hermitian matrix and z a complex 
Gaussian vector of n dimensions with mean m and covari- 
ance R. The matrix R is assumed to be non-singular, as 
it is a covariance matrix, is therefore positive definite. 

Expression ([6]) only makes sense if the determinant 
1 1 — juRQ| is non-zero. Proof of this fact is easy to find if 
the determinant is expressed as Yii=i (1 — 3 v &i) UM > where 
{Si} represents the eigenvalues of RQ, which must be real 
so that the determinant will be non-zero. These eigenval- 
ues, by definition, satisfy |RQ — 6{l\ = 0, which can be 
rearranged as (R^QR 1 / 2 - 6il\ = 0, where R 1 / 2 is the 
matrix square root of R, which is also a Hermitian positive 
definite matrix. Hence, as matrix R 1 / 2 QR 1 / 2 is Hermi- 
tian, the eigenvalues {Si} are real and therefore |I — jwRQ| 
is non-zero. 

The characteristic function <f> in ([6]) can be rearranged as 
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This expression, particularized to 2 x 2 quadratic forms, is 
applicable to the variable dk of ([2]) to obtain Other- 
wise, it is easy to show that for 2 x 2 non-singular matrices 
M] and M 2 , whose sum is also non-singular, the following 
property is fulfilled 
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Using this property in expression ([7]) particularized to vari- 
able dk, and remembering that Q is non-singular, gives 
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or, if S\ and 82 are the eigenvalues of matrix RQ, 
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It is easy to show that the eigenvalues of matrix RQ can 
be defined as 



tr(RQ) ± Vtr(RQ) 2 - 4 |RQ| 



(11) 



where tr(-) represents the trace. The positive sign of (|11[) 
has been chosen for eigenvalue Si and the negative sign 
for eigenvalue 82- In order to obtain non-trivial results for 
Pr{D < 0}, as mentioned in the previous section, matrix Q 
must be indefinite. Hence, as matrix R is positive definite, 
matrix RQ is indefinite. So, Si takes a strictly positive 
value and S2 a strictly negative value. 

Comparing (ITU1) and J3]) it is possible to find the following 
alternative definitions for parameters v±, v 2 , aifc, 0:2k, 



vi = l/Si 
v 2 = -I/S2 
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These definitions are now in a compact matrix format 
and second they make it easier to find little mistakes in 
original expressions. This mistakes are clarified in table HI 
where the original parameters are listed in one column and 
the parameters with mistakes are corrected in the other 
column. It is important to highlight that the detected mis- 
takes only affect parameters w and aik insofar as they are 
related to constant C . Thus, the original expressions are 
only wrong when the constant C is complex. 

New expressions of the probability Pr{D < 0} can be 
used taking into account (|12j) . These expressions only 
needs the two eigenvalues, 81 and S 2 , and parameters a 
and b, and they are summarized in table ITT] 

In most cases where the original expressions have been 
used up until now, the constant C was real, so the results 
obtained were valid. For example, in chapter 12 of [2], in 
the context of multichannel digital communication with bi- 
nary signaling, two types of processing at the receiver are 
considered, coherent and non-coherent detection. For the 
coherent detector the constant C is equal to 1/2 and for the 
non-coherent detector the constant C is equal to 0. Thus, 
in both cases the constant C is real and the mistakes in the 
expressions do not affect to the results. Another example 
is presented in [B], where Gray-code 16-QAM constellation 
is considered. In this case, the decision variables for the 
in-phase and quadrature components are defined with con- 
stant C equal to 1/2 and —j/2, respectively. Despite the 
fact that constant C is complex for the quadrature compo- 
nent, the symmetry of the problem was taken into account 
in [6j, making it possible to perform the calculations in 
another way in which C was real. Therefore, the results 
obtained in jBJ are valid but would be wrong if the problem 
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TABLE I 

Expressions of the error rate for binary signals 
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TABLE II 

Alternative expressions of the error rate for binary signals 
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were solved using the quadrature component variable deci- 
sion (C complex). A simple example where C is complex 
is shown in the next section. Again, the example serves 
to show the mistakes in the original expressions, revealing 
how they lead to incorrect results and showing how the new 
expressions yield correct results. 

IV. Example 

In order to justify the implications of the mistakes in the 
original expressions, a simple example where the parameter 
C is complex is presented in this section. The example 
chosen has the following decision variable 

D =exp(j^)* 1 y 1 * + exp(-^)A 1 *y 1 

r tt t 4 (is) 

=2Re{exp(j-)A 1 r l *}, 

so, L = 1 and the constants are identified as A = 0, 
B = 0, and C = exp(jj). The complex Gaussian variables 
X\ and Y\ are chosen to be independent, with mean and 
variance X\ = exp(jj), Y\ — 1, and E jj A"i — A\ | 2 j = 

E< 111 — Yi| } = 1, respectively. 



TABLE III 

Parameters calculation 



Parameter 


Proakis 


Corrected 


w 








Vl 


1 


1 




1 


1 


an 


2 


2 




2 





a 





1 


b 


V2 


1 


Pr{D < 0} 


0.18394 


1/2 



Taking into account the definitions of table HI all param- 
eters are calculated in both ways, with original definitions 
and with corrected definitions, as presented in table Hm 
Note that the value of Pr{D < 0} calculated with the cor- 
rected expressions is exactly 1/2, for which eq.(4.53)] 
has been used. As shown in table Mil the two types of def- 
initions give different values for the last three parameters 
and the final probability. The example is chosen so that 
the final probability can be found easily in another way as 
follows. 

The probability Pr{D < 0} can be directly calculated 
from the definition (|T3"]) of D, as 

Pr{£> < 0} =Pr{2|Ai| |Yj.| cos ( J + LX X - ZYi) < o} 

= Pr{sin(ZYi - Z) < 0} 
= Pr{sin([/) < 0}, 

(14) 



0.12 




Fig. 1 

Histogram of decision variable D 



Any non-zero mean complex Gaussian variable V = 
pexp(j#) with mean V = rexp(j'O) and standard devi- 
ation a has the following joint probability density function 

(pdf)mi 

(15) 

It is evident from (jT5j) that the marginal pdf pe(0) has even 
symmetry when its mean O is zero. 

Note that X\ and Y\ are non-zero mean complex Gaus- 
sian variables. Therefore, the variables £Y\ and Z are sta- 
tistically distributed as 0, and as they have zero-mean they 
have an even pdf. Moreover, the difference U between the 
independent variables ZYi and Z also has even symmetry. 
This implies that the variable sin(f) has an even pdf, be- 
cause it is an odd function applied to a variable with even 
pdf. Consequently, the probability Pi{D < 0} is exactly 
1/2. 

Figure Q] shows the histogram of the decision variable D 
obtained by simulations. Note the even symmetry of the 
histogram, confirming that Pi{D < 0} = 1/2, the same 
value obtained with the corrected expressions. 

V. Conclusions 

In p] is presented an analysis of complex Gaussian 
quadratic forms which gives some powerful expressions to 
study many digital communications systems in terms of 
the error probability. This analysis has been revised in 
this technical report, obtaining alternative expressions in 
a compact matrix format. Little mistakes in the original 
expressions have been detected which must be corrected 
to avoid errors in practical applications. Finally, a simple 
example has been presented which reveals the influence of 
these mistakes and confirms the validity of the new expres- 
sions. 
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